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Abstract: Oculopharyngeal muscular dystrophy (OPMD) is caused by trinucleotide repeat expansion mutations in Poly(A) 
binding protein 1 (PABPNl). PABPNl is a regulator of mRNA stability and is ubiquitously expressed. Here we investigated 
how symptoms in OPMD initiate only at midlife and why a subset of skeletal muscles is predominantly affected. Genome- 
wide RNA expression profiles from Vastus lateralis muscles human carriers of expanded-PABPNl at pre-symptomatic and 
symptomatic stages were compared with healthy controls. Major expression changes were found to be associated with 
age rather than with expression of expanded-PABPNl, instead transcriptomes of OPMD and elderly muscles were 
significantly similar (P<0.05). Using k-means clustering we identified age-dependent trends in both OPMD and controls, but 
trends were often accelerated in OPMD. We report an age-regulated decline in PABPNl levels in Vastus lateralis muscles 
from the fifth decade. In concurrence with severe muscle degeneration in OPMD, the decline in PABPNl accelerated in 
OPMD and was specific to skeletal muscles. Reduced PABPNl levels (30% to 60%) in muscle cells induced myogenic defects 
and morphological signatures of cellular aging in proportion to PABPNl expression levels. We suggest that PABPNl levels 
regulate muscle cell aging and OPMD represents an accelerated muscle aging disorder. 



INTRODUCTION 

Protein aggregation is a pathological hallmark in a large 
spectrum of late-onset neurodegenerative disorders. 
Single amino acid repeat expansions are among the 
common genetic cause of those disorders [1]. These 
disorders are recognized by the presence of salt- 
resistant insoluble bodies, resulting from accumulation 
of the expanded mutant proteins. Accumulation of 
aggregated proteins is associated with impaired protein 
homeostasis and the ubiquitin proteasome system 



(UPS), which advances protein aggregation [2]. The 
presence of disease-associated protein aggregation is 
often an indication of failure of cellular processes 
maintaining protein homeostasis [3]. Nevertheless, 
despite their evident occurrence in late-onset 
neurodegenerative disorders, a pathogenic role of these 
aggregates remains controversial [4]. 

Most of the expansion mutations create elongated 
polyglutamine (poly-Q) tracts, but a subset of diseases 
is caused by polyalanine tract expansions [5]. 
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Autosomal dominant oculopharyngeal muscular 
dystrophy (OPMD) is caused by a poly-alanine repeat 
expansion mutation in the gene encoding for Poly (A) 
Binding Protein Nuclear 1 (PABPNl) [6]. OPMD is 
characterized by progressive ptosis, dysphagia, and 
proximal limb muscle weakness that typically appear 
from the fifth decade [5]. Accumulation of expanded 
(exp)PABPNl into insoluble intranuclear inclusions in 
skeletal muscles is the pathological hallmark of OPMD 
[7, 8]. Because of the many commonalities in the 
protein aggregation process, it is often considered a 
paradigm for the neurodegenerative protein aggregation 
disorders with the unique advantage of having access to 
pre-symptomatic and symptomatic tissue. Naturally 
occurring inclusions of wild-type (WT) PABPNl have 
also been reported [9]. In contrast to aggregates of 
expPABPNl, those of the WT protein are not associated 
with disease. In muscle cell models both WT and 
expPABPNl form nuclear inclusions, but expPABPNl 
is more prone to aggregation [10]. This increased 
aggregation of expPABPNl can be explained, in part, 
by its reduced poly-ubiquitination and slower protein 
turnover compared with WT PABPNl [10]. 

The impact of aggregated PABPNl on cellular defects 
is not fully understood. In high overexpression 
situations expPABPNl, but not WT PABNPl, induces 
cell death [11, 12]. Therefore, a toxic gain-of-function 
was suggested for mutant PABPNl (reviewed in: [5]). 
In cellular models with low overexpression levels that 
are similar to levels of endogenous PABPNl, and 
therefore mimic better the situation in heterozygous 
patients, apoptosis was not observed [10, 13]. In 
myotube cultures with low PABPNl overexpression, 
deregulation of the ubiquitin proteasome system (UPS) 
was the most affected cellular pathway [10]. 
Deregulation of the UPS was also found to be most 
prominent in a cross-species study for OPMD [14], and 
reduced proteasome activity was found in the A 17.1 
mouse model for OPMD [15]. This suggests that 
disruption of cellular processes regulating protein 
homeostasis, rather than cell death, could cause muscle 
degeneration in OPMD. 

High overexpression of expPABPNl or depletion of 
PABPNl by RNA interference in muscle causes 
myogenic defects [16-20]. These models have been 
instrumental to demonstrate that manipulations in 
PABPNl expression levels cause muscle symptoms. As 
yet, however, it is unclear how ubiquitously expressed 
PABPNl causes predominant skeletal muscle pathology 
and why symptoms in OPMD initiate only after midlife. 
Here we investigated the hypothesis that OPMD 
represents accelerated aging myopathy. We compared 
genome-wide expression profiles from Vastus lateralis 



(VL) muscles of OPMD patients and elderly and 
identified significant similarities. In concurrence with 
severe muscle weakness in OPMD, age-dependent 
molecular changes are often augmented in OPMD 
muscles. We also identified an age-dependent decline in 
PABPNl levels. This decrease in PABPNl levels starts 
only from the fifth decade and is accelerated in OPMD. 
In culture, PABPNl down-regulation induces muscle 
cell aging and myogenesis defects that are proportional 
to a decrease in PABPNl levels. We suggest that 
PABPNl regulates muscle cell aging. 

RESULTS 

Transcriptome changes in OPMD are associated 
with initiation of symptoms 

Symptoms in OPMD often initiate with eye lid muscle 
weakness (ptosis) or and swallowing difficulties 
(dysphagia) [7], but muscle weakness in upper limbs was 
also reported [21]. In our patient cohort limb muscle 
weakness was often highly abundant (Table 1). Muscle 
strength in 33 genetically confirmed OPMD patients 
(Supplementary Table 1) was determined with the 
Medical Research Council (MRC) scale, where MRC <5 
represents clinical symptoms. An association between 
muscle strength and age was statistically assessed from 
1 1 muscle groups, from which significant association was 
found in 8 muscle groups (Table 1). Using the same 
muscle strength measures, in carriers of expPABPNl 
under 42 years-old muscle weakness was not observed 
(Supplementary Table 1). This suggests that age is major 
confounder in muscle weakness in OPMD. 

In OPMD VL muscles we identified genome-wide 
changes in gene expression [14]. Here we assessed 
whether these changes in OPMD are associated with 
with the expression of expPABPNl and/or with age. 
We first compared the transcriptome of OPMD patients 
with that of expPABPNl carriers at a pre-symptomatic 
stage. Expression profiles for each group were 
generated after subtraction to an age-matched control 
group: for the pre-symptomatic group mean age was 36 
±3.5 N=6 and 9 (carriers and controls, respectively) 
and in the symptomatic group mean age was 56 ± 5 
years N=9 and 13 (carriers and controls, respectively) 
(Supplementary Table 1). The genes found to be 
affected in OPMD (N=1572 genes; p-value < 0.01) 
were clustered into 73 functional Kyoto Encyclopaedia 
of Genes and Genomes (KEGG) pathways and were 
highly similar to those identified in OPMD animal 
models (Supplementary Table 2). With the same p- 
value cut-off only 473 genes were identified in the pre- 
symptomatic group and those did not enrich in any 
KEGG pathway. This suggests that gene deregulation in 
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pre-symptomatic have little functional impact. After 
applying the false discovery rate (FDR)-correction, 149 
genes were found to be deregulated in the symptomatic 
group while in the pre-symptomatic group no gene 
passed those stringent criteria. This analysis thus shows 
that in OPMD muscle transcriptome changes are 
associated with symptoms, while changes caused by the 
presence of expPABPNl are only minor with limited 
functional impact. Although more confidence of gene 
deregulation is obtained after FDR correction, since our 
dataset is relatively small we selected the deregulated 
genes with /7-value 0.01 for further analysis. 



Table 1. Muscle weakness in OPMD is age- 
dependent 



Muscle (group) 


Beta ofMRC 

age-associated 

1 


p-value 


Iliopsoas 


-0.07 (0.02) 


3.8E-04 


Biceps 


-0.03 (0.01) 


2.0E-03 


Hamstrings 


-0.04 (0.01) 


5.0E-03 


Feet 


-0.02 (0.01) 


5.0E-03 


Neck extensors 


-0.03 (0.01) 


6.0E-03 


Med. Gluteal 


-0.05 (0.02) 


2.0E-02 


Quadriceps 


-0.02 (0.01) 


3.0E-02 


Triceps 


-0.01 (0.01) 


4.8E-02 


Deltoid 


-0.03 (0.01) 


8.0E-02 


Hip extensors 


-0.08 (0.02) 


2.7E-01 


Neck flexors 


-0.01 (0.02) 


4.5E-01 



Regression analysis of MRC scores dependent on age (in 
years) of patients at examination. Betas (standard errors) 
are provided. P-values are adjusted for gender. N= 33; 
IVIean age in years: 60.4 (± 11.3); Females 49%. 



Expression profiles in OPMD muscles and in 
muscles of elderly are highly similar 

We first investigated the possibility of molecular 
similarities between OPMD muscle and muscle aging 
using functional groups. The affected KEGG pathways 
in OPMD were highly similar to those found in an 
independent study of normal muscle aging [22]. From 
the 75 KEGG pathways that are significantly affected in 
this muscle aging study only 3 (96% similarity) were 
not identified in OPMD (Supplementary Table 2). 

If molecular changes in OPMD are similar to these in 
muscle aging similarity in dysregulated genes and 
expression trends is expected. We compared OPMD 
expression profiles (mean age 56 ± 5) with those of 



elderly (>85 years) (Supplementary Table 1). It is 
relevant to note that OPMD patients were under 70 
years-old and therefore may not considered as elderly. 
VL muscles from healthy controls were sampled same 
as the OPMD biopsies and expression profiles were 
generated on the same platform. Expression profiles in 
OPMD and elderly were generated with the same 
control group (mean age 56 ± 5). Close to one third of 
the genes affected in elderly overlapped with those 
affected in OPMD (Figure lA), and the same 
directionality in fold-change was found for 84% of 
these genes. The majority of the overlapping genes 
(77%)) were down regulated (Figure lA). This 
observation is consistent with other studies that found a 
prominent transcriptional down-regulation in neurons 
from the elderly [23]. Gene Ontology (GO) analysis of 
OPMD- or elderly- regulated genes also revealed high 
similarities between significantly affected cellular and 
molecular functions (Table 2). 

We then investigated whether aging- or OPMD- 
regulated genes are correspondingly enriched in 
recognized aging concepts. We employed the literature 
association method [24] with the concepts: aging, 
muscle contraction, oxidative phosphorylation, insulin 
signalling, tumour growth factor (TGPP) signalling and 
the ubiquitin proteasome system (UPS). The literature 
association of every gene with each of the 6 concepts 
was determined and normalized association-weights 
were assigned. Association enrichment for the 
significantly regulated genes in each concept was 
evaluated with the Cumulative Distribution Function 
(CDF). The distribution of association-weights for the 
aging concept had a maximum score indicating that the 
common genes between OPMD and elderly indeed 
represent well-known aging genes (Figure IB). As 
expected, significant scores were also found for 
oxidative phosphorylation and muscle contraction. Next 
to the aging concept, the UPS also rated with a 
maximum score, suggesting that UPS genes have 
highest impact for both OPMD and elderly. The high 
significant score indicates that the association is not 
random {P = 4.3x10"^^ 8.1x10"^^ and 2.4x10"^^ for 
UPS, oxidative phosphorylation and muscle contraction 
respectively). In contrast, the distributions of weight- 
associations for genes in the insulin or TGF-P signalling 
pathways were not strong (Figure IB), and P-values did 
not differ from a theoretical random distribution. This 
analysis suggests that similar genes involved in the 
UPS, mitochondria and muscle contraction affect 
muscle weakness in both elderly and in OPMD. The 
genes affecting the insulin or TGF-P signalling 
pathways, however, differ between OPMD and elderly. 
To visualize the molecular similarity between OPMD 
and elderly we plotted the mean fold change for well- 
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Table 2. Top 10 affected GO-terms in OPMD and elderly 



OPMD 



GO term 



P-value Nr 



Elderly 



GO Term 



P-value 



Nr 



000568 l~spliceosome 

000558 l^collagen 

0044420'-^extracellular 
matrix part 

00065 11 ^ubiquitin- 
dependent protein 
catabolic process 

000628 l^DNA repair 

O000123'-histone 
acetyltransferase 
complex 

O016570'-histone 
modification 

00069 15'-apoptosis 

0006405~RNA export 
from nucleus 
O044429'-mitochondrial 
part 



I.IE- 
07 

0.002 



0.003 

0.004 

0.004 

0.007 

0.019 
0.026 
0.031 
0.044 



29 



31 

35 

10 

17 
59 
8 

54 



000568 1-spliceosome 

0044429-mitochondrial 
part 

000651 l^ubiquitin- 
dependent protein 
catabolic process 

0043292^contractile fiber 

0015629'-actin 
cytoskeleton 

00043 1 1-famesyltranstran 
s-ferase activity 

0033279^ribosomal 
subunit 

O006006'-glucose 
metabolic process 
0003678-'DNA helicase 
activity 

0043065^positive 
regulation of apoptosis 



5.9E-05 
l.OE-04 

0.006 

0.008 

0.009 

0.023 

0.027 
0.032 
0.040 
0.051 



22 
61 

27 

16 
28 
3 

15 
17 
7 

37 



Enriched GO terms were identified with the 
The nunnbers (Nr) of genes nnapped to each 



list of OPMD- or Elderly-affected genes (P-value <0.01). 
GO-ternn are indicated. 



OPMD 



Elderly 




C 2.5 



same direction 

down (77.5%) 
up (22.5%) 



■ OPMD 
□ Elderly 



Insulin signalling 

TGpp signalling 



Muscle contraction 



Oxidative 
phosphorylation 
Aging 
V 




0.2 0.4 0.6 0.8 

Normalized Association Weights 



3130 



ititit iO 



Figure 1. Similar expression profiles in OPMD and 
elderly muscles. (A) Venn diagrann shows the nunnber of 
affected genes (p <0.01) in OPMD, elderly overlapping 
genes. The percentage of sinnilar gene direction and the 
direction (up or down) are listed. For each group the nnean 
age ± standard deviation are denoted. (B) Cunnulative 
distribution function (CDF) plots show the distribution of 
normalized literature association-weights of commonly 
deregulated genes between OPMD and elderly with the 
concepts: Aging, Muscle contraction. Oxidative phospho- 
rylation, insulin signalling, TGFp signalling, and the 
ubiquitin-proteasome system (UPS). Arrowheads indicate 
the maximum association weights. The further the curve 
shifts to the right, the higher the associations of the 
affected genes with the indicated concepts are. (C) Bar- 
chart shows fold-change of selected known age-regulated 
genes in elderly and in OPMD. P-values are indicated: * P 
<0.01, ** P<0.005, *** P<0.0005. 
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known aging-regulated genes [25] (Figure IC). 
Together our analysis revealed a significant 
consistencyin molecular signatures that characterize VL 
muscles in OPMD patients and in elderly. 

Transcriptome of an OPMD mouse model clusters 
with muscle aging transcriptome 

We further studied the association of the OPMD 
transcriptome with aging using the transcriptome of the 
A17.1 mouse at 6-weeks old, as muscle weakness is not 
obvious at that age (Trollet et al, 2010). In the A17.1 
mouse expPABPNl is overexpressed under a muscle- 
specific promoter resulting in severe muscle weakness 
(Davies et al., 2005) and muscle atrophy at an early age 
[15]. A literature-aided association study of the 
differentially expressed genes in the A17.1 muscles of 6 
weeks-old mice showed that a large subset of these 
genes strongly associates with the Siging concept (Figure 
2A). The fold-change for the genes with a high 
association-weight score was large (Figure 2A), 
indicating that well-known aging genes are highly 
dysregulated in this mouse model. 



For an unbiased approach we performed an 
unsupervised meta-analysis with 104 microarray studies 
of muscle development and muscle disorders [26], 
including the transcriptome from muscle of 6-weeks old 
A 17.1 mice. With this analysis similar transcriptional 
changes were grouped into three main clusters (Figure 
2B). The transcriptome of the 6 week-old A17.1 mouse 
strictly grouped together with those related to skeletal 
muscle aging [22, 27], but not with datasets from other 
muscular dystrophies or myopathies (Figure 2B). This 
suggests that the high overexpression of expPABPNl in 
the OPMD mouse model could induce similar 
molecular changes as in normal muscle aging. 

Expression trends are accelerated in OPMD 

Close to 30 years separate the OPMD patients in our 
study from the elderly, but the remarkable overlap in 
expression profiles suggests that changes in OPMD 
initiate earlier or progress faster compared with normal 
aging. To identify expression trends k-means clustering 
was applied on the overlapping probes between OPMD 
and elderly. We compared trends per probe in a dataset 
derived from the healthy cohort age 17-89 years to 
those derived from expPABPNl carriers age 31-74, and 
subsequently, unique genes were identified using Entrez 
ID. Two clusters of significant age-regulated genes 
showing accelerated (N=35) or decelerated (N=39) 
trends in OPMD were identified (Figure 3 A, 
Supplementary Table 3). Examples of trends for 
selected genes are presented in Figure 3B. Among those 
are well-known regulators of aging such as the cell 
cycle regulators CDKNIA (p21) and OSBP, regulators 
of muscle cells and energy metabolism like LMODl and 
CHRNAl or FATl and PRODH. These genes with 
accelerated expression changes in OPMD could 
contribute to disease progression and muscle weakness. 

PABPNl levels significantly decreases with age and 
in OPMD 

The mRNA processing and spliceosome GO-terms were 
the most significantly affected in both OPMD and 
elderly datasets (Table 2). Noticeably, the majority of 
the dysregulated genes in this category were down 
regulated (Figure 4A). Network analysis suggests that 
central hubs are affected in both OPMD and elderly, 
including PABPNl (Figure 4B). 

Figure 2. The transcriptome of the OPIVID mouse model is highly associated with aging. (A) Volcano plot shows the distribution of 
significantly deregulated genes (P = 0.05; indicated with a dashed line) in 6 week-old A17.1 nnice against fold-change. Snnaller P-value and higher 
fold-change suggests a higher innpact in OPMD. For every gene a literature-associated weight with the 'Aging' concept is assigned, and a 
nornnalized association-weight is presented with a circle on a linear scale between 0.05 and 1, where 1 equals the highest association. (B) 
Hierarchical clustering arrangennents of 104 datasets in a literature-aided nneta-analysis. Shades of blue indicate degree of sinnilarities: fronn weak 
(white) to strong (dark blue). Three skeletal nnuscle aging-related datasets cluster together with the OPMD dataset of 6 week-old nnice 
(highlighted in red). The clusters associated with nnuscular dystrophies and other nnyopathies are highlighted in green and blue, respectively. 
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Since mutations in PABPNl are the genetic cause for 
OPMD, we next validated the decHne in PABPNl 
mRNA levels using RT-qPCR using an extended cohort 
of skeletal muscle biopsies. A significant and 
pronounced decrease in expression was found in OPMD 
compared with age-matching controls, while at the pre- 
symptomatic stage PABPNl levels did not significantly 
change (Figure 5 A). This suggests that symptoms in 
OPMD patients are associated with a decrease in 
PABPNl expression. 

Next, levels of PABPNl mRNA during aging were 
determined in VL from 74 healthy controls age 17-89 
using RT-qPCR. Fold-change was calculated from 
GUSB housekeeping gene and young control group at 
age 17-22 years. A significant decrease in PABPNl 
expression with age was observed with a quadratic 
model (Figure 5C). Using two sequential linear models 
the decline in this dataset was found from 43 years 
onwards, whereas between 17-42 years the change in 
expression levels was insignificant (Figure 5C). 
Together, our analysis reveals that PABPNl expression 
declines with aging from the fourth decade onwards and 
that this decline is accelerated in OPMD. 

PABPNl is expressed in every cell whilst symptoms in 
OPMD are predominantly restricted to a subset of 
skeletal muscles. To investigate whether the decrease in 
PABPNl expression can explain the tissue specificity in 
OPMD, PABPNl levels in blood samples of OPMD 
patients and controls was determined using RT-qPCR. 
PABPNl levels in blood were unchanged between 
OPMD patients and age-matching controls (Figure 5B), 
indicating that the decrease in PABPNl expression in 
OPMD is more prominent in skeletal muscles. An age- 
associated decrease in PABPNl levels was also not 
found in cross-sectional microarray studies from blood, 
parotid glands, kidney cortex and medulla (Table 3). In 
frontal brain cortex we identified a slight decrease in 
PABPNl expression (Table 3), which started only after 
the age of 70 years (Table 3). Also in Rectus abdominis 
muscles, where muscle aging is not so prominent [28], 
PABPNl expression did not change with age (Table 3). 
Together, this suggests that the age-associated decrease 
in PABPNl expression marks aging of human skeletal 
muscles that are typically affected in OPMD. 



Figure 3. Age regulated gene expression trends change 
faster in OPIVID. (A) Expression trends in controls (grey) and 
expPABPNl carriers (red). Sinnilar expression trends were 
identified witli l<-nneans clustering using probe ID. (B) Exannples 
of expression trends of 8 genes fronn clusters in A, in healthy 
controls (grey) or in exPABPNl carriers at pre-synnptonnatic and 
synnptonnatic stages (red). Pre-synnptonnatic and synnptonnatic 
stages of OPMD are indicated. 



PABPNl down-regulation (-DR) induces myogenic 
defects and muscle cell aging 

Knockdown of PABPNl in mouse muscle cells causes 
myogenic defects [18]. Since in vivo we found a 
decrease in PABPNl expression with age in both 
normal aging and in OPMD, we studied whether at with 
a moderate down-regulation of PABPNl myogenic 
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defects can also be observed. Three PABPNl shRNA 
clones were selected for functional studies in immor- 
talized human myoblast cultures using the lentivirus 
expression system for efficient delivery. Stable cultures 
were generated and down-regulation was evaluated with 
RT-qPCR and western blot analyses. Compared with 
controls (HI empty vector and non- transduced cells), the 



three PABPNl shRNA clones, 121, 122 and 123, led to a 
70%, 40% and 20% decrease in PABPNl expression, 
respectively (Figure 6Ai). This decrease was also reflect- 
ed in the residual protein levels (Figure 6Aii). Immuno- 
fluorescence microscopy with antibody to PABPNl 
revealed a decrease in the accumulation of nuclear 
PABPNl in the shl 21 -transduced cells (Figure 6Aiii). 




Figure 4 -Deregulation of the Spliceosome in both OPMD and elderly. (A) Bar-chart shows fold-change of spliceosome genes 
in OPMD and elderly fronn the transcriptonne studies. P-values are indicated: * P<0.05, ** P<0.005, *** P<0.0005. In bold are 
gene-hubs in the OPMD-affected spliceosonne network (B). A schennatic gene network presentation of OPMD-regulated genes 
that are grouped in the spliceosonne GO category. PABPNl is highlighted 
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Figure 5. RT-qPCR analysis of PABPNl expres- 
sion trends in OPMD and during muscle aging. 

(A) Box plot shows P>ABPA/lexpression in Vastus 
lateralis from expPABPNl carriers at a pre- 
synnptonnatic (pre-symp) or symptomatic (OPMD) 
stages, and age-matching control groups. (B) Box plot 
showsP>ABPA/lexpression in blood form OPMD patients 
and controls. (C) Scatter plot shows P>ABPA/lexpression 
in 78 healthy controls age 17-89 years. Male and 
female samples are indicated in black and grey, 
respectively. A quadratic fit is shown with a black line 

(a) , and linear fits are for the age groups: 17-42 years 

(b) or 43-89 years (c) are denoted in grey. The table 
summarizes p-values and Beta ± standard errors, 
which were calculated after gender correction. Fold- 
changes were calculated after normalization of GUSB 
housekeeping gene and to control group age 17-22 
years (A, C) or age matching controls (B). N denotes 
the number of samples in each group. 
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Table 3. Changes in PABPNl expression depends on clironological age are 
muscle specific 



Tissue 


Age 




Beta 


P-value 


Vastus lateralis 


17-42 y(N = 


41) 


-0.006 (0.009) 


0.37 


muscles 


43-89 y (N = 


34) 


-0.029 (0.006) 


<0.0001 


Frontal Brain Cortex 


26-69 y (N = 


17) 


0.002 (0.007) 


0.73 




70-95 y (N = 


13) 


-0.018 (0.008) 


0.04 


Blood 


42- 102y(N = 


= 150) 


0.001 (0.003) 


0.69 


Kidney Cortex 


27 - 92 y (N = 


-72) 


-0.001 (0.002) 


0.76 








-0.001 (0.002) 


0.42 








-0.003 (0.002) 


0.15 


Kidney Medulla 


29 - 92 y (N = 


= 61) 


-0.003 (0.002) 


0.11 








0.001 (0.002) 


0.76 








-0.004 (0.002) 


0.06 


Rectus Abdominis 


24-83y(N = 


81) 


-0.000 (0.003) 


0.94 








0.010 (0.007) 


0.13 








0.001 (0.003) 


0.64 


Parotid glands 


19-71y(N = 


= 13) 


0.000 (0.003) 


0.93 








0.003 (0.005) 


0.64 








-0.001 (0.005) 


0.86 



Betas (standard errors) and p-value are provided per probe. Values for three 
independent PABPNl probes are shown for datasets fronn Kidney cortex, Kidney 
nnedulla, Rectus Abdominis and Parotid glands. P-values are adjusted for gender. 
Significant changes are highlighted in bold. N indicates number of sannples. Age is 
indicates in years (y). 



Next, we investigated whether PABPNl -DR 
specifically affect expression levels of myogenic genes. 
We selected MYHl, DMD and CAV3 for RT-qPCR 
analysis, as their expression levels significantly 
decreased in both OPMD and elderly. The strongest 
decrease in expression of all three genes was found in 
the shl21 cultures with highest PABPNl -DR (Figure 
6B). Only a moderate decrease across all three genes 
was found in the shl22 culture, and expression levels 
did not differ from controls in the shl23 culture (Figure 
6B). The decrease in MHYl on mRNA level was 
confirmed on the protein level (Figure 6A). This suggests 
that levels of PABPNl -DR can induce a decrease in 



cell fusion. An apparent decrease in cell fusion was 
found in the strong PABPNl -DR, shl21, cultures 
(Figure 6Ci). We then determined rate of cell fusion 
across PABPNl -DR fused cultures index using a semi- 
automated and robust imaging procedure, where the 
percentage of nuclei within cells expressing MHCl was 
quantified. The rate of cell fusion was not affected in 
the mild PABPNl down-regulated culture, shl23, but in 
shl22 or shl21 cultures cell fusion significantly 
reduced (Figure 6Cii). Together the results here reveal 
that a decrease in PABPNl expression induces 
myogenic defects that are proportional to levels of 
PABPN 1 -down-regulation. 
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Figure 6. PABPNl-DR in human myotubes causes myogenic defects. Human myotubes were transduced with shRNA 
specific to PABPNl (slil21, slil22, or slil23) or HI ennpty vector. Non-transduced (NT) cells were used as controls. (A) i Bar- 
chart shows PABPNl mRNA expression in stably-transduced nnyoblasts. Fold change was nornnalized to GapDH housekeeping 
gene and to a non-transduced culture. Averages are of 6 biological replicates. 11 Western blot analysis of PABPNl, MHCl and 
nnuscle actin (MSA), as a loading control in shl21, shl22 or HI nnyotube cultures, ill Innnnunofluorescence of PABPNl (labelled 
with Alexa-594) and MHCl (labelled with Alexa-488) in shl21 or HI nnyotube cultures. Scale bar 10 |Linn. (B) Bar chart shows 
fold change of MHCl, DMD, and CAV3 in 121-, 122-, 123-, and HI- nnyoblast cultures. Fold change was nornnalized to GapDH 
and to a non-transduced culture. Averages are of 3 biological replicates. Significant down-regulation (P<0.05) is indicated with 
asterisks. (C) 1- innages of HI controls and the PABPNl down regulation shl21 fused cultures. Scale bar is 50 |inn. ii- Chart bar 
shows the percentage of nuclei in fused nnyotubes that express MHCl (cell fusion) in HI shl23 shl22 or shl21 nnyotube 
cultures. Averages and SD are from six replicates and the number of nuclei that were quantified per sample is indicated 
within each bar. Significant effect in PABPNl-DR cultures from control cultures (P<0.05 or P<0.005) is indicated with one or 
two asterisks, respectively. 



Aging cells undergo substantial cellular changes [29]. To 
investigate whether a decrease in PABPNl expression 
induces cellular characteristics of aging cells, we applied 
four different cellular assays. Reduced mitochondrial 
metabolic rate is often found in aged skeletal muscles 
[30]. Using JC-1 labeling, we found a decrease in 
mitochondrial metabolic rate in both myoblast and 
myotube cultures with PABPNl down-regulation 
compared with control cultures (Figure 7Ai). The effect 
of PABPNl-DR on the mitochondrial metabolic rate was 
more pronounced and more significant in myotube 
cultures compared with myoblasts cultures (Figure 7Ai). 
Importantly, a decrease in mitochondrial metabolic rate 
was restored in PABPNl-DR cultures after transduction 
with CFP-PABPNl lentivirus particles (Figure 7Ai). 
Noticeably, the decrease in metabolic rate, levels of 
myogenic genes and fusion index correlated with levels 
of PABPNl-DR. As levels of mitochondrial metabolic 
rate are decreased due to aging-regulated increase in 



oxidative stress [31], next we investigated whether 
PABPNl affect cell sensitivity to oxidative stress. 
Muscle cell cultures were treated with H2O2 and 
mitochondrial metabolic rate was determined after JC-1 
labeling. PABPNl-DR muscle cells are 1.5-fold more 
sensitive to H2O2 treatment compared with control cells 
(Figure 7Aii). In addition. Levels of PABPNl also 
affect rates of cell growth. PABPNl-DR cells have a 
slower growth rate compared with control, and restored 
after transduction with CFP-PABPNl lentivirus 
particles (Figure 7Bi). In PABPNl-DR cells 
heterochromatic foci, which mark senescent cell [32], 
were also detected (Figure 7Bii). Last, we also found 
that PABPNl-DR induces formation of fat droplets in 
muscle cell cultures (Figure 7C). Importantly, restoring 
PABPNl levels with CFP-PABPNl decreased the 
proportion of cells with fat droplets (Figure 7Ci). 
Together, these studies indicate that reduced levels of 
PABPNl induce muscle cell aging. 
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Figure 7. PABPNl-DR induces cellular aging in muscle cells. (A) PABPNl-down regulation induces a decrease in 
mitocliondrial nnetabolic rate, (i) Bar cliart sliows tine ratio of J-aggregates to monomers in JC-l-labeled myoblasts or 
myotube cultures. Cell cultures are of control (7304) or PABPNl-DR (shl22 and shl21) cultures, before (mock) or after 
transducion with CFP-PABPNl lentivirus particles. Averages and SD are from six replicates, and the number of cells per 
sample is indicated within each bar. Significant effect in PABPNl-DR from control cultures (P<0.05 or P<0.005) is 
indicated with one or two asterisks, respectively, (ii) Plot shows the ratio of J-aggregates to monomers in JC-l-labeled 
cultures after 1 hour H202, in concentrations as indicated in the chart. Averages and SD are from three replicates. (B) A 
decrease in cell growth is caused by PABPNl-down regulation, (i) Bar chart shows cell growth (24 hours) in control 
(7304) and PABPNl-DR (shl21) cultures, before (mock) or after transduction with CFP-PABPNl lentivirus particles. 
Averages and SD are from three replicates, (ii) Heterochromatic foci are formed in PABPNl-DR cells. Heterochromatic 
foci are visualized with DAP! (indicated with arrowheads), cells are visualized with Desmin. Scale bar is 20 |Lim. 



DISCUSSION 

Protein aggregation is a hallmark of a spectrum of late 
onset neurodegenerative disorders. A pathogenic or 
toxic role for aggregated proteins in general, and hence 
the functional relevance of expPABPNl aggregates, to 



disease pathogenesis is under debate [4, 33]. Since 
PABPNl aggregates deplete levels of the soluble 
protein [10], a pathogenic model that postulates the 
reduced bioavailability of PABPNl as a consequence of 
PABPNl aggregation, seems equally plausible. In 
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support, in vitro, the function of soluble PABPNl is 
concentration dependent [34], and in cell models, 
similar changes in PAS usage are induced by either high 
overexpression of expPABPNl (thus forming PABPNl 
aggregates and hence depleting the cell from soluble 
PABPNl) or PABPNl down-regulation [35]. In muscle 
cell models with similar overexpression levels of WT or 
expPABPNl the changes in PAS usage were limited 
arguing against a gain-of- function mechanism [36]. 
This suggests that the levels of soluble PABPNl rather 
than a gain-of function of mutant PABPNl affects 
muscle cell function. Here we found that the expression 
of PABPNl is dependent on age. Here we found a 
natural decrease in PABPNl mRNA levels in VL 
muscles starting from the fifth decade onwards and 
progresses more rapid in OPMD patients compared with 
controls. Depletion of PABPNl from muscle cells 
causes myogenesis defects [18]. Here we show that in 
addition to myogenesis defects, PABPNl down- 
regulation induces phenotypes of muscle cell aging, 
including reduced cell growth rate, heterochromatic foci 
formation, reduced mitochondrial metabolic rate, fat 
droplets formation and increased sensitivity to oxidative 
stress. In addition, we show that levels of myogenesis 
and mitochondrial metabolic rate correspond to levels 
of PABPNl expression. Here we provide the first 
indications that an age-regulated decline in PABPNl 
levels, rather than gain-of function of mutant PABPNl, 
may cause an increase in muscle weakness. However, it 
is remained to be determined whether a decrease in 
PABPNl levels with age is caused by increase in 
atrophy, increase in central nuclei, and/or increase in 
fibrosis, fat and inflammation in theses muscles. Further 
cell-cased assays should be implied to determine the 
change in PABPNl levels in different cell types. In 
OPMD, severe muscle weakness would be the result of 
depletion of soluble PABPNl due to protein 
aggregation and an aging-regulated decrease in 
expression levels. 

PABPNl is a regulator of mRNA stability [37]. 
Knockdown of PABPNl expression in cells causes a 
decrease in poly(A) tail length [18] and a switch from 
distal to proximal polyadenylation site (PAS) usage at 
3'-UTR [35, 36]. Thus, changes in PABPNl levels 
could directly affect genome-wide expression profiles. 
In previous studies we reported that in OPMD and 
models with expPABPNl overexpression 
transcriptomic changes are genome-wide (Anvar et al., 
2011). The expression profiles in muscles from OPMD 
and elderly are highly similar, despite the average age 
difference of 30 years, suggesting that common 
molecular mechanisms regulate these genome-wide 
transcriptional changes. The mRNA processing and the 
spliceosome machinery were significantly deregulated 



in both OPMD and elderly, implying an aging- 
associated impairment of mRNA processing in skeletal 
muscles. Alternative splicing has been previously 
suggested to contribute to transcriptional changes 
during cell senescence, due to changes in expression 
levels of splicing factors [38]. Alternative splicing is 
also thought to contribute to late-onset neuro- 
degenerative disorders [23]. The impact of mRNA 
processing to genome-wide expression levels during 
aging and the molecular details will in the near future 
likely be uncovered by next generation sequencing 
studies. 

PABPNl is ubiquitously expressed but symptoms in 
OPMD are predominately restricted to a subset of 
skeletal muscles. We found that in OPMD patients 
PABPNl expression decreases in affected muscles but 
not in blood. It is important to note that expression 
profiles generated from muscle biopsied reflects cell 
heterogeneity in the actual biopsy, and those expression 
profiles may not be specific for muscle cells. However, 
age-regulated changes in PABPNl levels were not 
found in six different tissues. In brain cortex a slight, 
but signiflcant decline in PABPNl expression was 
found from the eight decade. In VL muscles the decline 
in PABPNl starts from the fourth decade. This suggests 
that a decrease in PABPNl levels is associated with 
aging-associated muscle weakness and that this process 
is accelerated in OPMD as judged from the trend that 
the levels of a substantial number of genes changes 
faster in OPMD muscles compared to muscles of 
healthy controls. Together this suggests that OPMD 
should be regarded as a disorder of accelerated muscle 
aging. 

METHODS 

Pre-symptomatic carriers, patients and controls and 
clinical investigation. We invited 33 patients from a 
cohort of 10 genetically conflrmed Dutch OPMD 
families. All symptomatic and pre-symptomatic carriers 
were confirmed by genetic analysis. A standard 
neurological examination including MRC score was 
performed as described previously in part of this cohort 
[39]. Muscle biopsies were obtained from 22 patients 
and 74 healthy controls (Supplementary Table 1). All 
biopsies are from VL muscles and were collected using 
the Bergstrom needle procedure. The biopsies froze 
immediately in liquid nitrogen. Blood was taken from 
all patients for genetic confirmation. Asymptomatic 
descendants of these patients were invited to participate 
without disclosure of their individual histological and 
genetic results. Eighteen of them underwent a muscle 
biopsy and genetic testing. Twelve descendants who did 
not carry the mutation were classified as controls. 
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whereas 6 descendants who did carry the mutation 
where identified as pre-symptomatic. This study was 
approved by the local Ethical Committee (CMO nr. 
2005/189) and written informed consent was obtained 
from all patients. Biopsies from healthy controls were 
collected at Canisius-Wilhelmina Hospital, Nijmegen, 
The Netherlands. 

Microarray and Statistical Analyses 

Datasets: The human and mouse OPMD microarray 
datasets are publicly available at GEO repository under 
the accession numbers GSE26605 and GSE26604, 
respectively and have been previously published [14, 
15]. In all datasets genome-wide expression profiles of 
skeletal muscles from OPMD were compared to 
controls. 

Data Processing: Quantile normalization was applied 
on the microarray raw dataset and data quality was 
assessed by the principal component analysis. 
Differentially expressed genes between two age-groups 
were identified by applying hierarchical linear model, 
using limma package in R[40] at a P-value <0.01. 
OPMD and pre-OPMD regulated genes were identified 
from age-matching control groups; and elderly 
regulated genes were identified from controls age 56± 5 
years. Probe annotation was carried out using 
illuminaHumanvSBeadID (human) and illumina- 
MousevlBeadlD (mouse) R packages. 

The meta-analysis in mouse was carried out on 104 
microarray datasets from various organisms as 
described in [26]. The 103 datasets are listed in [26], 
including the A17. 1 6 week dataset. 

The association-weights between genes and biological 
concepts were mined using Anni 2.1 [24]. To compare 
between gene groups, the association weights were 
normalized to the maximum association weight. Genes 
with literature association <0.1 and P>0.05 (-logio>1.3) 
in elderly and OPMD were excluded. Cumulative 
Distribution Function (CDF) analysis was used to 
examine the distribution of association weights, and the 
Kolmogorov-Smimov (KS) test was used to identify 
distributions that significantly differ from a theoretical 
distribution, threshold of P <10"^. The CDF of Genei is 
defined as the proportion of genes with association 
weight less than or equal to that of Genei, 

A k-means clustering of affected probes was used to 
identify age-regulated expression trends. However, in 
order to optimize the clustering arrangements, average 
Silhouette (S^yg) values are calculated for each cluster in 
Matlab. Clustering arrangement of partitions with 



Ss^yg<0.6 were reiterated until the criteria has met. 
Maximum number of clusters was set to 20 to avoid 
overly complex clustering arrangement due to the size 
of the set. The cluster centroids were used to provide 
summarized age-dependent expression patterns for each 
cluster. 

Mapping of affected genes to functional groups (KEGG 
or GO) was carried out was carried out with the 
globaltest [41] or with DAVID functional annotation 
clustering tool[42]. Gene network analysis was 
performed with STRING [43]. 

Statistical analyses were carried out with the SPSS 
software (IBM) and Matlab. 

Cell culture and Lentivirus transduction. The human 
7304 immortalized myoblasts [44] were propagated in a 
medium containing DMEM+20% Fetal Calf Serum 
supplemented with an equal volume Skeletal Muscle 
Cell Media (PromoCell, Heidelberg, Germany) at 37 °C 
under 5% CO2. Cell fusion was carried out in a medium 
containing DMEM+5% Horse Serum. 

The shRNA in lentivirus expression vectors 121 
(TRCN0000000121), 122 (TRCN0000000122) and 
(TRCN0000000123) 123 were obtained from Sigma- 
Aldrich. An empty vector, HI, was used as a negative 
control. Lentivirus production and transduction was 
previously as described[45]. Cells were cultured with 
viruses (MOI -'25) overnight, followed by medium 
refreshing. Transduced cells were maintained in the 
presence of 5jLig/ml puromycin, which was omitted for 
fusion experiments. PABPNl down-regulation was 
determined using RT-qPCR. Six independent 
transduction experiments were performed. Cell fusion 
was carried out in triplicates and cell fusion index was 
determined by dividing the number of nuclei in 
myotubes to the total number of myotubes. 

The lentiviruses, CFP-PABPNl and H4-CFP, were 
previously described in ([13, 45], respectively). 
Transduction of each virus was performed in control 
(7304), shl21 or shl22 cultures. Only cultures with 
transduction efficiency >70% were included in the 
studies presented here. 

RNA Analysis. All human biopsies were stored at -80 
before RNA extraction. RNA extraction and RT-qPCR 
were performed as described in [14]. RNA extraction 
from culture cells was performed with NucleoSpin® 
RNA II, MACHEREY-NAGEL, according to the 
manufacturer protocol and RNA was directly subjected 
for RT-qPCR using Cyber Green, Invitrogen. 
Expression levels were calculated according to the 
AACT method, and were first normalized to GUSB 
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housekeeping gene for human samples and GAPDH for 
cell culture. Fold change was calculated from controls. 
The statistical significance was determined with the 
Student's t-test. The list of primers used in this study is 
provided in Supplementary Table 4. 

Western blot, immunofluorescence, JC-1 and Nile-red 
staining and image quantification 

Western blots: Total proteins were extracted from fused 
cells and subjected to western blot analysis as described 
in [10]. Detection of the first antibodies was conducted 
with the Odyssey Infrared Imaging System (LI-COR 
Biosciences) and suitable secondary antibodies. 

Immunofluorescence: of fused cells was carried out on 
cells seeded on plastics or on collagen-coated glass 
plates as described in [10]. Cell cultures were fused for 
7 days and after fixation and permeabilization were 
incubated with first and subsequently with the 
appropriate secondary-conjugated antibodies. The 
Alexa 488-, Alexa 430- or Alexa 594- conjugated 
secondary antibodies against primary antibodies were 
obtained from Molecular Probes (Invitrogen). 

Primary antibodies used in this study are: mouse 
monoclonal anti-muscle actin (MSA) (Novocastra), 
(1:2000) and anti-Myosin MF20 (1:500) (Sigma- 
Aldrich). PABPNl was recognized with the 3F5 llama 
single chain [46] antibody and rabbit-anti-VHH. 

Nile-red staining: Nile red is a fiuorescent stain for 
intracellular lipid droplets [47]. Cells were seeded on 
glass and after 24 hours were incubated with 1 |liM Nile- 
Red (Invitrogen) for 20 minutes. Cultures were washed 
and mounted with citifiour containing DAPI (500 
ng/ml). Nile-Red was imaged with a Cy5 filter. 

JC-1 staining: JC-1 was used to measure the 
mitochondrial membrane potential [48]. 1 |liM JC-1 
(Invitrogen) and Hoechst (Invitrogen) for nuclear 
counter staining in growth media and directly added 
into myoblasts or 7 days fused cultures for 20 minutes. 
Cultures were washed and living cells were directly 
imaged with the Array Scan VTI HCA (Thermo 
Scientific) using XF53 filter set. Nuclei were segmented 
with the nuclear counterstain and the ratio of red to 
green signal was quantified from compartments 
surrounding each nucleus using the Compartmental 
Analysis V4 BioApplication (Thermo Scientific). Plots 
were generated in Excel. 

Cell fusion: Fusted cultures that were immunelabeled 
with MF20 and counterstained nuclei were scanned 
with the ArrayScan VTI HCA (Thermo Scientific) 



using UV to detect nuclei and FITC to detect the MF20 
signal. The percentage of nuclei within MF20 objects 
was identified with the Colocalization V4 
BioApplication (Thermo Scientific). Plots were 
generated in Excel. 
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